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Abstract 

Using the previous Belle measurement of the inclusive photon energy in B — > X s j decays, we 
determine the first and second moments of this spectrum for minimum photon energies in the 
B meson rest frame ranging from 1.8 to 2.3 GeV. Combining these measurements with recent Belle 
data on the lepton energy and hadronic mass moments in B — > X c tv decays, we perform fits to 
theoretical expressions derived in the IS and kinetic mass schemes and extract the magnitude of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V c b, the 6-quark mass and other non- 
perturbative parameters. In the IS scheme analysis we find \V c b\ = (41.56 ±0.68(fit) ±0.08(tb)) x 
1CT 3 and mf = (4.723 ± 0.055) GeV. In the kinetic scheme, we obtain | V cb \ = (41.58 ± 0.69(fit) ± 
0.08(r B ) ± 0.58(th)) x 10~ 3 and mf n = (4.543 ± 0.075) GeV. 

PACS numbers: 12.15.Ff,12.15.Hh,12.39.Hg,13.20.He 
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I. INTRODUCTION 



The most precise determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ment | V c b | [1] are obtained using combined fits to inclusive B decay distributions [2-5]. These 
analyses are based on calculations of the semileptonic decay rate and spectral moments in 
B — > X c lv and B — > X s ^ decays in the frameworks of the Operator Product Expansion 
(OPE) and the Heavy Quark Effective Theory (HQET) [2, 6-8], which predict these quan- 
tities in terms of \V c b\ and a number of non-perturbative heavy quark (HQ) parameters 
including the 6-quark mass m&. 

Analyses combining measurements from different experiments [2, 3] quote the most pre- 
cise numbers for \V c b\ and m&. However, as the correlated systematic uncertainties are not 
precisely known, there is some concern that uncertainties are underestimated. In this anal- 
ysis, we have chosen the opposite approach and perform fits to the data from the Belle 
experiment only. In addition, we use two independent sets of theoretical expressions, de- 
rived in the IS [2] and kinetic mass [7, 8] schemes respectively, to test the compatibility of 
these two frameworks. 

The present document is organized as follows: Sect. II describes the measurement of 
the first and second moment of the inclusive photon energy spectrum in B — > X s -y, (E^) 
and ((Erf — (Erf)) 2 ), using the Belle measurement of this decay in Ref. [9]. In the previously 
published analysis the first and second moments were obtained for one value of the minimum 
energy threshold, namely E min = 1.8 GeV. Here we report additional measurements with 
-^min = 1.9,2.0,2.1,2.2,2.3 GeV, and perform a re-evaluation of the systematic error. In 
Sect. Ill we use these data together with the recent Belle measurements of the lepton energy 
and hadronic mass moments in B — > X c iv decays [10, 11] to extract \V c b\ and m& using 
theoretical expressions derived in the IS and kinetic mass schemes. 

II. MOMENTS OF THE B X s ~f PHOTON ENERGY SPECTRUM 
A. Review of the Belle B — ► X s 7 Measurement 

The analysis described in Ref. [9] uses e + e~ — > T(4S) — > BB events equivalent to 140 fb^ 1 
of integrated luminosity (ON sample) and 15 fb _1 taken 60 MeV below the T(45 l ) resonance 
energy (OFF sample). Photon candidates with energy greater than 1.5 GeV as measured in 
the T(4S I ) rest frame are reconstructed. Vetoes are applied to photon candidates with high 
likelihood of originating from ir° or rj decays to two photons. 

In general, the background of photons from the e + e~ — > qq continuum is dominant. It is 
suppressed with event shape variables used as the inputs to two Fisher discriminants [12]. 
The first discriminant distinguishes spherically-shaped BB from jet-like continuum events 
and includes the Fox- Wolfram moments [13], the thrust calculated using all particles detected 
in the event including and excluding the candidate photon, and the angles of the correspond- 
ing thrust axes with respect to the beam and candidate photon directions, respectively. The 
second discriminant is designed to exploit the topology of B — > A s 7 events by utilizing the 
energy sum of detected particles, which is measured in three angular regions bounded by 
cones that are subtended from the direction of the candidate photon in the T(4S) frame; 
defined as 0° - 30° (forward), 30° - 140° (middle), and 140° - 180° (backward). 

After these selections are applied, the remaining continuum background is removed by 
subtracting scaled OFF data from the ON data set. Backgrounds in BB events, including 
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FIG. 1: Left: raw photon energy spectrum in the T(4S) frame; right: photon energy spectrum 
after background subtraction and efficiency correction where the inner error bars are the statis- 
tical uncertainties and the outer error bars show the total errors, which include the systematic 
uncertainties. These plots are reproduced from Ref. [9]. 

photons from 7r° and t] (veto leakage), other real photons (mainly from u, rf, and J/ip), 
clusters in the calorimeter not due to single photons (mainly electrons interacting with 
matter, K\ and ft) and beam background, are estimated from Monte Carlo (MC) simulation 
(Fig. 1). 



B. Moment Measurements 

We calculate the truncated first and second moments, (E^) and ((i? 7 — (_E 7 )) 2 ), of the 
efficiency corrected spectrum in Fig. 1 for minimum photon energies ranging from 1.8 to 
2.3 GeV. The following corrections are applied to these moments: The non-zero B meson 
momentum in the T(4S*) rest frame changes the first moment of the photon energy by 0.2% 
and adds a Doppler broadening of 0.006 GeV 2 to the second moment; the finite energy 
resolution, uncorrected in Fig. 1, causes a broadening of the spectrum and increases the 
second moment by 0.004 GeV 2 ; the 100 MeV binning in Fig. 1 increases the second moment 
by 0.0008 GeV 2 . 

The above corrections assume a symmetric photon energy distribution, and do not ac- 
count for expected and known asymmetries in the true spectrum and detector response, 
respectively. To account for these effects an additional bias correction, derived from a MC 
simulation, is implemented. The B — > model contains decays of the form B — > K*~f, 
where K* is any known spin-1 resonance with strangeness 5 = 1. The relative amounts of 
these decays are adjusted by matching the total photon spectrum to the theoretical model 
of Ref. [14]. The bias correction, calculated as the difference of the true moment and the 
moment measured in the B —>■ X s ^ MC simulation once all aforementioned corrections are 
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TABLE I: Residual bias correction to (E^) and ((E^ — (E^)) 2 ) as a function of E, 



E min (GeV) 


A(E 7 ) 


A<(£ 7 -<£ 7 )) 2 > 


1.8 


+2.0% 


0.0% 


1.9 


+1.6% 


-0.4% 


2.0 


+1.2% 


-7.1% 


2.1 


+0.8% 


-17.4% 


2.2 


+0.2% 


-35.3% 


2.3 


-0.3% 


-57.9% 



applied, is listed in Table I. 

C. Systematic Uncertainties 

The error bars of the efficiency corrected spectrum of Fig. 1 show the total error includ- 
ing the systematic uncertainty related to the scaling of the MC background samples (sizable 
in the first energy bins). In the calculation of the moments, we consider also the follow- 
ing sources of systematic uncertainty: uncertainty in the OFF data scaling factor; possible 
difference in ON and OFF data selection efficiencies; uncertainty in the BB data/MC cor- 
rection; we vary by ±20% the background from 77', u and bremsstrahlung; uncertainty of 
the 77 veto efficiency; we consider an alternate signal MC that favors high-mass resonances 
decaying into high-multiplicity final states, where the fraction of ^Kir final states, some- 
what overestimated in our default sample, matches our previous measurement [15]; and we 
vary the photon detection efficiency in both signal and background samples by its measured 
uncertainty (±2.3%). 

We also assign systematic uncertainties to the corrections applied to the moments: an 
alternate energy resolution correction that neglects the lower energy tail in the resolution is 
implemented and the difference is assigned as systematic uncertainty; a ±100% uncertainty 
on the binning correction for the second moment is assigned; we also implement a ±50% 
variation on the bias correction for the first moment while for the second moment the 
correction is re-calculated using the alternate signal MC sample. 

The total systematic uncertainty on each moment measurement is obtained by summing 
the aforementioned contributions in quadrature (Tables II and III). 

D. Results 

The measurements of the first and second moments of the photon energy spectrum in 
B — > Xg'-f for minimum photon energies ranging from 1.8 GeV to 2.3 GeV are shown in 
Table IV and Fig. 2. Our results agree with the data from CLEO [16] and BaBar [17]. 

The statistical and systematic errors on the first and second moments at E min = 1.8 GeV 
are slightly different from the values quoted in Ref. [9] and supersede our previously published 
values. The change in the uncertainties is due to the use of the toy MC approach and to 
the additional contribution from the uncertainty in the bias correction. 
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TABLE II: Systematic uncertainties contributing to the first moment (E~) as a function of the 
lower energy threshold E m i n in GeV. 



£min (GeV) 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


MC scaling 


0.021 


0.012 


0.006 


0.003 


0.002 


0.001 


OFF scaling 


0.004 


0.001 


0.001 


0.002 


0.002 


0.002 


ON/OFF efficiency 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


BB data/MC correction 


0.005 


0.003 


0.002 


0.001 


0.000 


0.000 


other 7s in BB 


0.010 


0.004 


0.002 


0.000 


0.000 


0.000 


77 veto efficiency 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


signal MC 


0.004 


0.004 


0.004 


0.004 


0.004 


0.003 


7 efficiency 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


bias correction 


0.022 


0.018 


0.014 


0.009 


0.002 


0.003 


total systematic 


0.033 


0.023 


0.016 


0.010 


0.005 


0.004 



TABLE III: Systematic uncertainties contributing to the second moment ((£7 — (Ey)) 2 ) as a 
function of the lower energy threshold E m i n in GeV 2 . 



E miQ (GeV) 




1.8 




1.9 




2.0 




2.1 




2.2 




2.3 


MC scaling 





0060 





0027 





0009 





0003 





0001 





0001 


OFF scaling 





0018 





0010 





0006 





0005 





0005 





0004 


ON/OFF efficiency 





0000 





0000 





0000 





0000 





0000 





0000 


BB data/MC correction 





0010 





0004 





0001 





0000 





0000 





0000 


other 7s in BB 





0024 





0008 





0002 





0000 





0000 





0000 


77 veto efficiency 





0003 





0001 





0000 





0000 





0000 





0000 


signal MC 





0007 





0005 





0004 





0003 





0002 





0000 


7 efficiency 





0003 





0001 





0000 





0000 





0000 





0000 


energy resolution 





0020 





0020 





0021 





0022 





0023 





0024 


binning 





0008 





0008 





0008 





0008 





0008 





0008 


bias correction 





0068 





0040 





0026 





0005 





0004 





0011 


total systematic 





0099 





0055 





0036 





0024 





0025 





0028 



The correlations between the different moment measurements are estimated using a toy 
MC approach: Starting from the efficiency corrected spectrum in Fig. 1 we create new 
spectra by generating values of a Gaussian random variable for the contents of each bin, 
where the mean and standard deviation of the Gaussian correspond to the bin yield and 
its uncertainty in the original spectrum. The moments and their fluctuations with respect 
to each other were measured for each generated spectrum, and finally averaged to yield the 
covariance matrix, from which the uncertainties due to statistics and systematics scaling 
were obtained. The covariance matrix was also obtained from systematic variations due to 
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TABLE IV: Measurements of (E~) and ((E~ — (E^)) 2 } as a function of the minimum photon energy 
E m in- The first error is statistical and the second is systematic. 



£ mi „ (GeV) 


<£ 7 ) (GeV) 


((E,-(E,)) 2 ) 


1.8 


2.292 ±0.027 ±0.033 


0.0305 ± 0.0079 ± 0.0099 


1.9 


2.309 ± 0.023 ± 0.023 


0.0217 ± 0.0060 ± 0.0055 


2.0 


2.324 ±0.019 ±0.016 


0.0179 ± 0.0050 ± 0.0036 


2.1 


2.346 ±0.017 ±0.010 


0.0140 ± 0.0046 ± 0.0024 


2.2 


2.386 ± 0.018 ± 0.005 


0.0091 ± 0.0045 ± 0.0025 


2.3 


2.439 ± 0.020 ± 0.004 


0.0036 ± 0.0045 ± 0.0028 



2.55 

2.5 
2.45 

2.4 
2.35 

2.3 
2.25 

2.2 









Stat 
Tota 


error 
error 
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FIG. 2: Measurements of (E~) and {(E^ — (E^)) 2 ) as a function of the minimum photon energy 
E min (AE 2 = ({EL, - (E 7 )) 2 )). 

the aforementioned corrections to the moments. The method assumes 100% correlation of 
any two truncated moments due to any single systematic variation. The covariance matrices 
that are derived from statistics and systematics are added to yield the overall covariance 
matrix, from which the correlations between any of the truncated moments are deduced. 
Tables V— VII show the correlation coefficients derived from this study. 

To cross-check these moment measurements, we extract the moments from the Kagan- 
Neubert (KN) photon spectrum [14] tuned to fit our data [18] (m 6 (KN) = 4.62 GeV, 
fi 2 (KN) = 0.40 GeV 2 ). We generate the photon spectrum in the rest frame of the B meson 
with these parameters and extract the moments in the range E m i a = 1.8, . . . , 2.3 GeV. The 
results are plotted in Fig. 3 along with the moment measurements presented here. We find 
very good agreement between these independent methods. 
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TABLE V: Correlation coefficients between the (E~) measurements. 



E ■ 

(GeV) 


1.8 


1.9 


<£ 7 
2.0 


) 

2.1 


2.2 


2.3 


1.8 


1.00 


0.79 


0.68 


0.56 


0.38 


0.22 


1.9 




1.00 


0.82 


0.70 


0.52 


0.33 


(E y ) 2.0 






1.00 


0.86 


0.67 


0.47 


2.1 








1.00 


0.84 


0.65 


2.2 










1.00 


0.86 


2.3 












1.00 



TABLE VI: Correlation coefficients between the (E^) and ((E^ — (E^)) 2 } measurements. 



E m {n 




((Ej- 


(E,)) 2 } 






(GeV) 


1.8 


1.9 2.0 


2.1 


2.2 


2.3 


1.8 


-0.46 


-0.18 -0.01 


0.04 


0.01 


-0.01 


1.9 


-0.06 


-0.21 0.05 


0.12 


0.10 


0.07 


(J5 7 ) 2.0 


-0.14 


0.15 0.12 


0.23 


0.20 


0.17 


2.1 


0.27 


0.37 0.43 


0.42 


0.39 


0.34 


2.2 


0.38 


0.55 0.67 


0.75 


0.66 


0.61 


2.3 


0.43 


0.63 0.79 


0.91 


0.88 


0.79 



TABLE VII: Correlation coefficients between the ((E^ — (E^)) 2 ) measurements. 



Emm 






((E 7 - 


(E,)) 2 ) 






(GeV) 




1.8 


1.9 2.0 


2.1 


2.2 


2.3 




1.8 


1.00 


0.72 0.63 


0.49 


0.39 


0.30 




1.9 




1.00 0.83 


0.71 


0.61 


0.52 


m- 


2.0 




1.00 


0.89 


0.80 


0.71 


-(E 7 )) 2 ) 


2.1 






1.00 


0.96 


0.91 




2.2 








1.00 


0.97 




2.3 










1.00 
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FIG. 3: Cross-check of the moment measurements. The moment measurements presented here 
are compared to the moments predicted in the Kagan-Neubert prescription, tuned to fit our data 
(AE* = ((E, - <£ 7 » 2 ». 

III. EXTRACTION OF \V cb \ AND m b FROM INCLUSIVE B DECAYS 
A. Experimental Inputs 

Belle has measured the partial branching fractions A3 and the first, second, third and 
fourth moments of the truncated electron energy spectrum in B — > X c tv, (Ei), ((Eg—(Eg)) 2 ), 
{{En — (Ee)) 3 ) and ((Eg — (Eg)) 4 ), for nine different electron energy thresholds (E min = 0.4, 
0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8 and 2.0 GeV) [10]. This analysis uses T(4S) -> BB events 
equivalent to 140 fb" 1 of integrated luminosity. The hadronic decay of one B meson is fully 
reconstructed and B — > X c ev decays of the other B are selected by requiring an identified 
electron amongst the particles remaining in the event. 

In addition, Belle has measured the first, second central and second non-central moments 
of the hadron invariant mass squared (M|) spectrum in B -> X c tv, (M x ), ((M% - (M|)) 2 ) 
and (Mx), for seven different lepton energy thresholds (E min = 0.7, 0.9, 1.1, 1.3, 1.5, 1.7 and 
1.9 GeV) [11]. This analysis is also based on 140 fb _1 of T(4S) data. Again, one B meson 
is fully reconstructed and a charged lepton (electron or muon) from the decay of the other 
B is required. The hadronic X c system is reconstructed by summing the 4-momenta of the 
particles remaining in the event. 

The measurements of the first and second moments of the photon energy spectrum in 
B — > Xg'-f, (Ery) and ((E 7 — (E y )) 2 ), have been described previously in this document. They 
are available for six different photon energy thresholds (E min = 1.8, 1.9, 2.0, 2.1, 2.2 and 
2.3 GeV). 

Hence, there are a total of 71 Belle measurements of inclusive spectra in B decays avail- 
able for the global analysis [19]. The measurements actually used in the IS and kinetic mass 
scheme fit analyses are shown in Table VIII: We have excluded measurements that do not 
have corresponding theoretical predictions; measurements with high E m[n cut-offs (i.e., elec- 
tron energy and hadronic mass moments with E min > 1.5 GeV and photon energy moments 
with -E m i n > 2 GeV) are not used to determine the HQ parameters, as theoretical expressions 
are not considered reliable in this region [8, 20]; finally, we have also excluded measurements 



9 



TABLE VIII: Experimental inputs used in the IS and kinetic mass scheme analyses. Both analyses 
use a total of 25 measurements. 





Measurements used 


Lepton energy moments {E™) 


n = 
n = 1 
n = 2 
n = 3 


-Emin = 0.6, 1.0, 1.4 GeV 
£ min = 0.6, 0.8, 1.0, 1.2, 1.4 GeV 
£ min = 0.6, 1.0, 1.4 GeV 
£ min = 0.8, 1.0, 1.2 GeV 


Hadronic mass moments (Mj^ 1 ) 


n = 1 
n = 2 


£ min = 0.7, 1.1, 1.3, 1.5 GeV 
E min = 0.7, 0.9, 1.3 GeV 


Photon energy moments (E~) 


n = 1 
n = 2 


E min = 1.8, 2.0 GeV 
E min = 1.8, 2.0 GeV 



where correlations with neighboring points are too high as these measurements do not add 
new information to the fit and introduce numerical problems such as negative eigenvalues 
of the covariance matrix. 

The value of | V c b\ is dependent on the average lifetime tb of neutral and charged B mesons. 
In the following analyses we use tb = (1.585 ±0.006) ps based on Ref. [21] and assume equal 
production of charged and neutral B mesons. 

B. IS Mass Scheme Analysis 

1. Theoretical Input 

The parameters appearing in the OPE depend on the choice of the mass scheme, i.e., the 
definition of mj. The IS scheme eliminates the 6-quark pole mass by relating it to the mass 
of the T(15*). Truncated spectral moments in B — > X c lv have been derived in this scheme 
up to 0{l/m\) [2]. The theoretical expressions are of the form 

(X) Einin = X™ + X&A + X^A 2 + X^A 3 + X^Xi + X^AX 1 + X^X 2 + X®AX 2 + X® Pl 
+X^p 2 + X^n + X^r 2 + X^r 3 + X^t 4 + X^e + X^e 2 BLM + X^eA , 

(1) 

where (X) stands for any experimental observable in Table VIII and X^\ i = 1, . . . , 17, 
are perturbatively calculable coefficients that depend on E min . The computations include 
radiative contributions of 0(e) and 0(e BLM ), the so-called BLM contribution at 0(e 2 ). The 
HQ parameters are A at leading order, Ai and A2 at 0{l/m 2 ), and t±, t 2 , T3, T4, p\ and 
p 2 at 0{l/m\). The CKM magnitude |V^,| enters through the predictions of the partial 
semileptonic branching fractions, 

G 2 F m b 2 

A#£ min = 1Q27r 3 \Vcb\ VQEDT-B(X)AB,E min , (2) 

where m is the IS reference mass, m = mx(is)/2, G 2 F m 5 / (192tt 3 ) = 5.4 x 10~ n ps -1 , 
^qed = 1.007 and (X)AB,E niin is an expression of the form of Eq. 1. 
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The analysis in the IS mass scheme determines a total of seven parameters: \V c b\, A, Ai, 
t~i, 72, T3 and p\. Following the prescriptions in Ref. [2], T4 is set to zero and the measured 
B* — B and D* — D mass splittings allow us to constrain some of the HQET matrix elements 
in Eq. 1: A2 = 0.1227 — 0.0145Ai and p 2 = 0.1361 + r 2 . The parameter A is the difference 
between the fe-quark and the reference mass, A = mx(is)/2 — m\ . We will present our results 
in terms of m£ s . 



2. The Fit 

The expressions in the IS scheme are fitted to the data using the x 2 minimization tech- 
nique and the MINUIT program [22]. The covariance matrix used in the fit takes into 
account both experimental and theoretical uncertainties. Following the approach in Ref. [2] , 
an element of the combined experimental and theoretical error matrix is given by 

a ij = a i a 3 C ij > ( 3 ) 

where i and j denote the observables and Cjj is the experimental correlation matrix element. 
The total error on the observable % is defined as 



* = y(0 J + W»0 J + (Bi/2) J 
for the nth hadron moment , 

d = yji*?*) 2 + {Af n {m B /2Yf + {BJ2Y 
for the nth lepton moment , 

*i = \I{°D 2 + {Af n {m B /2YY + {B.J2Y 

for the nth photon moment , (4) 

and /o = fi — 1, fi — 1/4 and f% = l/(6\/3). Here, a^ xp are the experimental errors, 
are the coefficients of the last computed terms in the perturbation series (used 
to estimate the uncertainty on the uncalculated higher order perturbative terms), and A 
is a dimensionless parameter that contains different theoretical uncertainties (uncalculated 
power corrections, uncalculated effects of order (a s /47r)AQ CD /m^, and effects not included 
in the OPE, i.e., duality violation). For lepton and hadron moments, we fix A = 0.001 [2]. 
For photon moments, the factor A is 0.001 multiplied by the ratio of the energy difference 
from the endpoint, relative to that for E min = 1.8 GeV, to account for the increase in shape 
function effects as one limits the allowed region of the photon spectrum. 

As the fit does not provide strong constraints on the l/m\ parameters, we add the 
following extra terms to the \ 2 function, 

v 2 -1° \(0)\<ml, 

X (o) ~ 1 i^z i//n\i ^ ™3 W 



m\) /Ml |(C)|>m 



where (m x ,M x ) are both quantities of O(Aqcd), and (O) are the matrix elements of any 
of the 0{l/m\) operators in the fit. For the central value of the fit, we take M x = m x = 
500 MeV [2]. 
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TABLE IX: Result of fit in the IS mass scheme. The er(fit) error contains the experimental and 
theoretical uncertainties in the moments. The (j{tb) error on \V cb \ is due to the uncertainty in the 
average B meson lifetime. In the lower part of the table, the correlation matrix of the parameters 
is given. 





\v cb \ (10- 3 ) 


m b (GeV) 


Ai (GeV 2 ) 


Pi (GeV 3 ) 


n (GeV 3 ) 


r 2 (GeV 3 ) 


r 3 (GeV 3 ) 


value 


41.56 


4.723 


-0.303 


0.067 


0.125 


—0.101 


0.125 


<r(fit) 


0.68 


0.055 


0.046 


0.030 


0.005 


0.056 


0.005 


o"(tb) 


0.08 














\v cb \ 


1.000 


-0.121 


0.003 


0.195 


0.008 


-0.432 


0.021 


m b 




1.000 


0.893 


-0.137 


-0.002 


-0.509 


-0.006 


Ai 






1.000 


-0.410 


-0.041 


-0.429 


-0.045 


Pi 








1.000 


0.009 


-0.533 


0.028 


T\ 










1.000 


0.005 


0.000 


T2 












1.000 


0.007 


T"3 














1.000 



TABLE X: Stability of the fit in the IS mass scheme. The different setups are explained in the 
text. Setup (d) corresponds to the default fit. 



Setup 


X 2 /ndf. 


\V cb \ (10- 3 ) 


m b (GeV) 


Ai (GeV 2 ) 


(a) 


6.4/14 


41.55 ±0.80 


4.718 ±0.119 


-0.308 ± 0.092 


(b) 


5.6/18 


41.28 ±0.86 


4.699 ± 0.060 


-0.491 ± 0.084 


(c) 


16.6/18 


41.10 ±0.54 


4.666 ± 0.046 


-0.341 ± 0.031 


(d) 


7.3/18 


41.56 ±0.68 


4.723 ± 0.055 


-0.303 ± 0.046 



3. Results and Discussion 

The results for the fit parameters are given in Table IX. Using the measurement of the 
partial branching fraction at -Emm = 0.6 GeV, we obtain for the semileptonic branching 
fraction (over the full lepton energy range) Bxju — (10.60 ± 0.28)%. A comparison of the 
measured moments and the IS scheme predictions is shown in Figs. 4 and 5. 

We have verified the stability of the fit by considering the following variations (Table X): 
(a) by repeating the fit only for the B — > X c lv data (21 measurements); (b) by releasing the 
m x constraint on the higher order parameters; (c) by repeating the fit with all theoretical 
uncertainties set to zero. In Table X the default fit corresponds to setup (d). Figure 6 shows 
the Ax 2 = 1 contour plots for the fits corresponding to setups (a) and (d) in Table X. 
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FIG. 4: Comparison of the measured electron energy moments and the IS scheme predictions 
(upper row), and difference between the measurements and the predictions (lower row). The error 
bars show the experimental uncertainties. The error bands represent the theory error. Filled circles 
are data points used in the fit, and open circles are unused measurements. 

C. Kinetic Mass Scheme Analysis 

1. Theoretical Input 

Spectral moments of the lepton energy and hadronic mass in B — > X c lv decays have been 
derived in the kinetic mass scheme up to 0(l/m\) [7]. Compared to the original publication, 
the theoretical expressions in our fit contain an improved calculation of the perturbative 
corrections to the lepton energy moments [23] and account for the i? min dependence of the 
perturbative corrections to the hadronic mass moments [24]. For the photon energy moments 
in B — > Xgj, the (biased) OPE prediction and the bias correction have been calculated [8]. 

These expressions depend on the following set of non-perturbative parameters: the b- and 
c-quark masses m& and m c , /i 2 and \j? G at 0(1/ ml) and (? D and p\ s at 0(l/m\) [25]. In 
our analysis, we determine these six parameters together with the semileptonic branching 
fraction (over the full lepton energy range) Bx c tv The total number of parameters in the 
fit is thus seven. 

The CKM magnitude |V^| is calculated using the following expression [6], 

= ( T m - X£u)^^] V2 x [1 + 0.30(a s - 0.22)] 
0.0417 V °- 105 / 

x [1 - 0.66(m fe - 4.6 GeV) + 0.39(m c - 1.15 GeV) 
+0.013(/4 - 0.4 GeV 2 ) + 0.09(p^ - 0.1 GeV 3 ) 

+0.05(/4 - 0.35 GeV 2 ) - 0M(p 3 LS + 0.15 GeV 3 )] , (6) 
where F(B — > Xiv) is the semileptonic width of the B meson. 
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FIG. 5: Same as Fig. 4 for the measured hadronic mass and photon energy moments and the 
IS scheme predictions. 



CM 
> 

a -o.2 



-0.3 



-0.4 



l 1 1 i 1 1 i 1 1 

Belle 

- 1S scheme 


i i | i i i i 










i 





4.6 



4.7 



— Q 0.043 



0.042 



4.8 4.9 

m b (GeV) 



0.041 



0.04 



1 


Belle 

1S scheme 


- 0X c lv 

_ 0X c lv+X s y 





4.6 



4.7 



4.8 4.9 

m b (GeV) 



FIG. 6: Ax 2 = 1 contours for the fit to all moments and the fit to the B — > X c £i/ data only. 
The Fit 

As in the IS scheme case, the fit is performed using the x 2 minimization technique and the 
MINUIT program [22]. The covariance matrix used is the sum of matrices corresponding 
to the experimental and theoretical uncertainties. The theoretical covariance matrix is 
constructed following the recipe in Ref. [7]: 

The non-perturbative uncertainties (i.e., the uncertainties related to the l/m& expansion) 
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are evaluated by varying /z 2 and p G (p% and p\ s ) by ±20% (±30%) around their "nominal" 
values of p\ = 0.4 GeV 2 , p 3 D = 0.1 GeV 3 , p? G = 0.35 GeV 2 and p\ s = -0.15 GeV 3 . 
The perturbative uncertainties (i.e., the uncertainties related to the expansion in a s ) are 
estimated by varying a s within ±0.04 (±0.1) around the central value of 0.22 (0.3) for the 
lepton and photon energy (hadronic mass) moments. The difference in the treatment of a s 
for the hadronic mass moments is due to the fact that the calculation of the perturbative 
corrections to these moments is less complete. 

The theoretical uncertainty in the moment predictions is the quadratic sum of these 
different contributions. The theoretical covariance matrix is then constructed by treating 
these errors as fully correlated for a given moment with different E min while they are treated 
as uncorrelated between moments of different order. 

For the moments of the photon energy spectrum, we take 30% of the absolute value of the 
bias correction as its uncertainty. This additional theoretical error is considered uncorrelated 
for moments with different E m [ n and different order. 

The experimental data from B* — B mass splitting and heavy quark sum rules constrain 
the parameters p? G and p\ s to 0.35 ± 0.07 GeV 2 and —0.15 ± 0.1 GeV 3 , respectively. We 
account for these constraints by adding the following additional terms to the \ 2 function, 

(p 2 G - 0.35 GeV 2 ) 2 /(0.07 GeV 2 ) 2 + (pis + °- 15 GeV 3 ) 2 /(0.1 GeV 3 ) 2 . (7) 

To calculate \ V c b\ using Eq. 6 and properly account for the correlations of the HQ parameters, 
we make \V c b\ a free parameter of the fit, calculate T(B — > XJLv) with Eq. 6 and add the 
following term to the \ 2 function, 



T(B -> X c iv) 



TB) 2 /a 2 . (8) 



The uncertainty o accounts for the experimental uncertainty in tb and an additional 1.4% 
theoretical uncertainty in extracting \V c -b\ using Eq. 6 [6]. We have verified that this method 
of calculating \V c b\ does not change the fit result for the other parameters. 



3. Results and Discussion 

The results of the fit in the kinetic mass scheme are shown in Table XL The value of the 
X 2 function at the minimum is 4.7 for 25 — 7 degrees of freedom. The semileptonic branching 
fraction Bx c iv is found to be (10.49 ± 0.23)%. The comparison of the measurements and the 
predictions in the kinetic scheme is shown in Figs. 7 and 8. 

We have repeated the fit using the B — > X c iv moments only, excluding B — > X s ^ data 
(Table XII). Figure 9 shows the A% 2 = 1 contour plots for the fits corresponding to setups 
(a) and (b) in Table XII. 



IV. SUMMARY 

We have determined the first and second moments of the photon energy distribution in 
B — > X s ^ decays, (E^) and ((E y — (i? 7 )) 2 ), for minimum photon energies in the B meson 
rest frame ranging from 1.8 to 2.3 GeV using the measurement of this spectrum published 
in Ref. [9]. The results are given in Table IV. We have also evaluated the (statistical and 
systematic) self- and cross-correlations between these measurements (Tables V— VII). 
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TABLE XI: Result of fit in the kinetic mass scheme. The cr(fit) error contains the experimental 
and theoretical uncertainties in the moments. The p{tb) and cx(th) errors on \V cb \ are due to the 
uncertainty in the average B meson lifetime and the limited accuracy of Eq. 6, respectively. In the 
lower part of the table, the correlation matrix of the parameters is given. 





\v cb \ (io- 3 ) 


m b (GeV) 


m c (GeV) 


A (GeV 2 ) 


P% (GeV 3 ) 


A G (GeV 2 ) pis (GeV 3 ) 


value 


41.58 


4.543 


1.055 


0.539 


0.166 


0.362 


-0.153 


er(fit) 


0.69 


0.075 


0.118 


0.079 


0.040 


0.053 


0.096 




0.08 














o-(th) 


0.58 














\Vcb\ 


1.000 


-0.371 


-0.316 


0.511 


0.493 


-0.166 


0.073 


m b 




1.000 


0.988 


-0.783 


-0.702 


-0.178 


-0.187 


m c 






1.000 


-0.771 


-0.715 


-0.262 


-0.108 


A 








1.000 


0.777 


0.205 


0.080 


P% 










1.000 


0.108 


-0.158 


A G 












1.000 


-0.103 


Pls 














1.000 



TABLE XII: Stability of the fit in the kinetic mass scheme. Setup (a) uses the B — ► X c iu data 
only; setup (b) corresponds to the default fit. 



Setup 


X 2 /ndf. 


\v cb \ (io- 3 ) 


m b (GeV) 


A (GeV 2 ) 


(a) 


4.2/14 


41.51 ±0.99 


4.573 ±0.134 


0.523 ±0.106 


(b) 


4.7/18 


41.58 ±0.90 


4.543 ± 0.075 


0.539 ± 0.079 



In the second part of the present document, we have combined these measurements with 
recent Belle data on the lepton energy and hadronic mass moments in B — > X c £u decays [10, 
11] to extract \V cb \, m b and other non-perturbative parameters using theoretical expressions 
derived in the IS [2] and kinetic [7, 8] schemes. 

The fits give consistent values of \V cb \ in the two schemes. In the IS scheme analysis we 
find |K&| = (41.56 ± 0.68(fit) ± 0.08(r B )) x IO" 3 and mf = (4.723 ± 0.055) GeV. In the 
kinetic scheme, we obtain jv^l = (41.58 ± 0.69(fit) ± 0.08(r B ) ± 0.58(th)) x 10~ 3 and mf a = 
(4.543±0.075) GeV. Note that the m b values can only be compared after scheme translation. 
The fit results using only the B -> X c tv data are \V cb \ = (41.55 ±0.80(fit) ±0.08(r B )) x 10~ 3 
and mf = (4.718 ±0.119) GeV in the IS scheme, and \V cb \ = (41.51 ±0.80(fit) ±0.08(r B ) ± 
0.58(th)) x 10~ 3 and mf n = (4.573 ± 0.134) GeV in the kinetic scheme (see Tables X and 
XII). 

The CKM magnitude \V cb \ and the fe-quark masses m£ m,is have been extracted with 
values that are consistent with previous determinations [2-5]. In the IS scheme \V cb \ has 
been measured with 1.6% precision. This is the most precise determination by any single 
experiment so far [4, 5]. 
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FIG. 7: Comparison of the measured electron energy moments and the kinetic scheme predictions 
(upper row), and difference between the measurements and the predictions (lower row). The error 
bars show the experimental uncertainties. The error bands represent the theory error. Filled circles 
are data points used in the fit, and open circles are unused measurements. 
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